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ABSTRACT _ 


This report presents a detailed description of a computer 
program to calculate cumulative damage of notched structural members 
subjected to arbitrary spectra. The Sequence Accountable Fatigue 


| Analysis computer program develops its sequence sensitivity by = 1% 


tracking residual stresses local to a notch throughout the/spectrum 


of loads. Residual stress relaxation analysis is included to increase _ 


% 


the generality of the results. An example spectrum and resulting 


cumulative damage analysis are illustrated. 
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SYMBOLS” 


Residual stress 

minea: local stress level | 

Minimum local stress level... 

Yield stress 

Equilibrium component of the residual stress 
Local strain, total | | - 
Elastic component of total strain 


Plastic component of total strain 


Maviaus apples stress level 


Mini mum applied stress level 


Mean applied stress level, Saget Suzn)/2 


Elastic stress concentration Factor 


_ Damage 


Integer describing the level number 


Equilibrium period, number of cycles for the local 
stresses to return approximately to the equilibrium 
conditions following an overload | 

Residual stress relaxation constant 

Relaxation function exponents | 

Number of cycles | 

Modulus of elasticity 


Number of cycles of Life at a given stress or. | 
strain cycling level | 


strain intercept at one reversal on a log 


Slope of the log ¢,~ log life curve 


ce ~log life curve 
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INTRODUCTION 


Cumulative danage a ah based upon the local stress-strain iz: 
behavior at a notch appear to be reasonably successful in anticipating 
trends in fatigue life behavior of notched specinens subjected ta Hedi 
spectrum Loading (1- ~6),. ' The type. of behavior that usually occura te 
that peak tensile loads tend to increase the fatigue life od Peak: oe 
coi reasive: lads tend to decrease the life of notched structureg_ bax 
‘compared to structures experiencing s0nd spectra not navsne those 
peaks (GS, 6), Local b behavior analyses, Buck as those developed by 
Smith ”) and Neuber (8) | , help to espiain this phenomenon as being 
a result of the tensile peak load creating a | compressive recdaan 
“stress at the notch and, conversely, the compressive peat aK 
a tensile residual stréss. The change in ‘life occurs because ‘phe aon 


residual stress state modi fies the subsequent damage accumulation ae 


The Sequence Accountable Fatigue Analysis computer Program wag 
developed to incorporate the local stress-strain epproacht with a ee 
recent residual stress relaxation analysis “) in order to Amprave oo 
the sequence pensrenvary of cumulative damage analysis. Thig 7 | 
technical memorandum presents ane details of the resultant computey | “i 
‘program and an ‘example of its, use. The correlation of predictions sre 
with this analysis to actual results of tests experiencing spectryy : 


. 2 loading is presented = Potter @ . and Potter and Noble , : 
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“PROGRAM OUTLINE 


The Sequence Accountable Fatigue Anaiyeta traces the stress-strain 
behavior local to a notch throughout an applied load spectrum and 
. enuciiaeas the auuaee based on the jecai experience. The computer 
program is divided generally into the four parts or modules outlined 
in Fig. 1. 

The basic input data for the material, specimen gecucty: fatigue " 
- behavior auebteies and spectrun, are developed in Module I. Module I 
is discussed further in Section III. Module IT paias the input infor- o 
nation and determines the local stress-strain behavior. Module IIl | 
references the Range Pair Counting Method Subroutine to eyele « count ees 
the local stress spectrum. Module IV determines ‘the damage in oa | 
local stress-strain oacerins | | . 

The basic sie a used in Modules It, Til and 1v are presented | 
below. — | | ia | 

Module II ~ Local Stress-Strain Behavior - . The analysis Saad during 
the cere pe atom of the local stress behavior during the spectrum 
; of iiedaine is a conbination of analyses developed by Smith @, ce 
ae (8) and Potter (6) Smith's simple analysis indicated that | 
she residual stress cola be approximated by assuming that ene initial 
etress-strain- behavior was elastic upon unloading following plastic = 


flow. Thus, the pestaiat: stress could be Keercurated snows.” the i o 
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in Fig. 2. 


aupiaa eauereed yielding). 


| oe min; 
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maximum local stress and the maximum applied stress as in Eq. 1 and 


Cis = Cmax, - ReSmax, a: (1) 


The Sequence Accountable Fatigue Analysis computer program currently 


incorporates elastic-perfectly plastic stress-strain behavior, tere 


fore, Yi. is equal to" "Ene yield stress. ‘For the cycles. timediately 


7 following the peak stress, the. residual stress “determined ‘in Eas Le . 
nodi fies the elastic solution: as shown in Eqs. 2 and 3. (provided that 
‘the following maximum anpiied: stress is gest than Sina , and that phere. — 


is no change in. the residual stress due to a. minimum applied atrpss 


| Cmax ~ Cres; _1 ida, i (2) 
er 7 Smaxy * Sresi_) * ae. eae aoe _ oe sae re 
the analysis developed by Neuber 8) has been extended to eye}te 
loading by Wetzel (2) and Wetzel, Morrow and Topper ( ° and used by 


| 1,4-6 ee 
many others ~ primarily to determine local stress-strain 7 


behavior. It is used in this program only to calculate plastic scpetna 


| occurring when the residual stress undergoes a | step change. The one 


plastic s strain calculation routine is accessed ently when the hae oe 
} 


it, + 
terms in Eqs. 2 and 3 exceed tensile or compressive yiel@ 


stress teverss respectively. ‘Figure 3 illustrates ae, 
of the plastic strain. Brae 5 | 
#2). 


The local stress-strain sehuerse, according to Wetzel ‘is- aaa abe 


- related to the applied load by Eq. 4 


ore #, «, s_ ye ai, | = arc ae 


2 piastic Stren can be Zoune by seneeaeeias a elastic come 
ponent rrom the eee strain. 


fp 7 SE > Se” (Ke Smax) “/E-Smax - Onax/E 
Therefore, the ae strain associated with om is given in ma. 5. 
| | i | neg 
“Pp; = (KS max, ) 2/Eoys ~ oy /E . ; ea See. @). 

if a see eee stress existed DESer to. this plastic strain | 


excursion, the plastic strain associated with that prior excursion 


: i subtracted from Eq. 53 as shown in Eq. 6. 


EDs = (KeSnax,) 2 Nays = (ys ~ Sresy_ o/ */ 80 oe = 6) : 

A similar calculation is made for plastic strains occurring during. 

the minimum stress peak. " # 
eco (6) 


In the analysis developed by Potter the residual stress 


cyclically relaxes toward zero or an equilibrium residual stress as 


shown in Fig. 4 according to Eq. ; 


Sresyer, 9 = (Sresy. oT ocanpg) 20g 1800 1)) Baki 
cease 
‘The Noy term, the Equilibrium Period, is dependent upon the. 

plied stress and the Residual Seeds Relaxation Constant. 


NEP, (C1/(K, 8, aes RSpean *)). 2 2 (8) ; 


‘the Residual Stress Relaxation conatent:,, Cl, has not yet been 


-xperimentally penned but should bé a constant for a material. - 


Module EIT - Cycle Counting Method 


Arter the local stress and sianelc strain behavior is calculated, 


. .ttad stress spectrum is Range Pair Counted using a ‘computer Program 


(11) 


Module IV - Damage Calculation — 


_ Damage is calculated separately for ‘the plastic strain excursions 


and the elastic stress spectrum. The damage is determined from the 


, conventional D “DN calculation. Damage from each of the plastic preatn ‘ 


cycles is determined from the Coffin-Manson ql ne 


(Dy = 1. -/Ne, = i I epgleg? Ve | 


"Damage from the elastic stress ee is determined in a. similar manner. | 


ee S-N data in a Modi fied Gosduan Diagram format. ‘Damage is 
bet summed, and failure. ‘of the coupon is’ defined as the event occurring when 


the summed damage equals unity. 


| The maximum and minimum focal stress levels are , sequentially ‘compared to 


; Sa wee re oc Co 2 ee Hed any sate aha, ie 
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INPUT DATA REQUIREMENTS © °° 


In general, each spectrum analyzed will require slightly different — 
programming in order to get the load history into a usable format for 


the core program. The basic. program requires : a certain family oe in- 


formation barous any analytical predictions can be nade. Appendix Lt 


contains a program listing for the ‘Sequence. Accountable Fatigue Anslysis, 
The subroutine CORE which accesses the subroutines having to do with RPCM, 


the Range Pair Counting Method, contains the basic analysis. ‘Subroutine | 


SAL reads the data input and then. references subroutine CORE. ‘The * 
‘subroutine SAL shown is one in which a block of cyeles is repeated with 
epeloual cycles. A list of the ca data cards and the ‘resulting . : 
cqciyels is et ven in Appendix ak - | : 
The specific data ‘requirements are given below. 
2 Stress-Strain Behavior - The stress-strain behavior is presuned 


to be elastic - perfectly plestse with the tensile yield stress ‘being 


equal to eine conpesaive yield stress. The yield stress value used is 


wii average of the monotonic behavior generally being above the 0.2% on 7. 


Siete value and balow the engineering ultimate strength. 

2. Residual Stress Relaxation | - The residual. stress relaxation | 
ctehavior of Eq. 7 and 8 ia characterized by cl, the Residual Stress a 
aaink Constant and El and E2, the atasdeton “equation caenies 


ee Residual Stress Relaxation Constant, cl, has not yet been adequately 


an eemiied: It should be a material senpatey if the estayecloa: function - 


a 


ay ree 


eS - . 


is correct and must be assumed. A reasonably accurate estimate of the 


Residual Stress Relaxation Constant for eimai material falls in the : 


range of 5-20 x 10° (cycles) (xsi). Further apex imentarton on the 


| pare of the ae should develop a ci isebis ‘for his set of conditions . 
until actual measurement of residual stress relaxation behavior defines ~ 
‘the relaxation. function and constants. At present El and £2 are eel as - 


| “to be equal to 1. 0. 


= Specimen Geometry ~ ‘The elastic K. value df avai lable) is | 


“entered into the analysis. If that eats is not available ‘then an 


estimate from some other method may be ised In certain cases, a value may 
be determined from a constant amplitude fatigue | test of a similar struc- 
ture by fitting several values of Ke to the analysis and deterring | 


the best correlation as is done with the Ke solution. Once a qtrese 


concentration factor, Kes is determined for a specinen, that value is 


not ee feo test-to-test : ee the same coupon configuration, 


4. Load Multiplier - Different spectra are presented for analysis 


in. different manners. Some data are presented in percent of maxiqym 


stress, eehets in terms of nominal poe and others in terms of pending 


moment... The value of the load multiplier defines the nominal, L atrags 


history. 


Ae “Cumulative Damage Analysis - The eae from the range-patyed 


elastic stress spectrum is determined by calculating a simple it yalue _* 


é 


- for each level and | accumulating the total. The Ne value is determined 


Z "from unnotched coupon S-N data in the es Goodman leah cae 


eee Se oes oe ee i ay rae ae 


Tae program ‘requires the input of four second order equations describing — 
the maximum and niaimum stress levels at lives of 10", 10°, 10° and 107 
cycles. The coefficients of the equations are aes by ase square 
fitting the S-N data presented in the form of Eq. 9. - os . 

| ax MD) Sai ie +B(I)S - + C(I) a - | = (9) te 
as A, B, and C aoetttetents for several typical materials are presented 
in Appendix Iv. The S-N data sieee was derived from various sources. 


ag The C coefficients s correspond. to = 


— usually oe the MIL ~HDBK~ oA 
the maximum stress level at zero to ‘maximum applied stress conditions ae 
on the unnotched coupons. — | | | 

The damage from che piasete cea cycles is determined using the 
Coffin-Manson relation to calculate the Nee value. ihe ‘conventional 
plastic strain intercept at one reversal and the Ey - life slope values 
are used in the sealy ete: Specific measured vature from the literature 
| are, aaea when available and typical values es they are not available. 
6. Analysis or Test Spectrum | - The last information needed is ‘the order 


and magnitude of application of the spectrum used in the test.” 
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OUTPUT OPTIONS 


i ee hk el ne a oe 


‘The Computer Progran petate the following output in the. process: 


of ene analysis. 


1. Maximum 208 minimum applied stress “and local stress eeepenee: 


through the spectrum. Also printed ¢ out is the residual stress, equilib- = 


| rium stress, applied cycles; and the equilibrium period. - a 


2. The elastic local stress history as ; input into the Range F Pair. 


"Subroutine and the resulting Range Paired spectrum: 


3. The plastic strain occurrence during the spectrum and the damage 


associated with Bach strain reversal. 


y 


4 The accumulated damage associated with the plastic strains. . 


5. The Range Paired elastic stress spectrum and the damage 


, associated with each level. | 


6. The accumulated “damage auuoetated with the current block of | 


loading including the plastic seratn damage and the total cemage since the 


initiation of cycling. 


At the option of the analyst, che can seit out all the above itens 


or only two. The IPRINT value controls what data is printed. 


—I£. IPRINT a all six items are > printed for each flight or block. 


i eb ip che ee ASSEN, ped ime pani in, ee eee 


If IPRINT = 2,. all items except 2. above are Re Or tees 


If IPRINT 3, only items 4. and 6. above are printed. 


ee 
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Coefficient of the x tem. (See A.) 
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A stress used ‘in the calculation of plastic strain, 


The constant. (See AD 


Inverse of the Coffin-Manson Np ea 


an ee ee 


The number of cycles in an Anterval. 0 oe 


The calculated camber ‘of cycles expected to be indicated 
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The Residual Stress: Relaxation Constant (see, ENEP y 


7 Damage. 


Decimal or real value of intecsens K after conversion, a8 


A portion of a least-squares-nethod solution, : 


The difference between residual stress and’ enue tht 


residual stress. - (RES (I) - EQRES) 


A portion of | a Jeast-squares-method : solution, Pec 


A variable used in the calculation of the nunber f cyples” 
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The mmber of applied cycles at a load level. > 


‘The ratio of the number of applied cycles to the number — 


of cycles to failure. ( ENN/ CYCLES ) 
LCF strain intercept. 
Equilibrium residual stress. 


An exponential function depicting the elaarion "ge | 


| residual stress. 


The name ‘of a "routine calling for the exponential value 


of a eres 


' An exponent. The power of 10 which indicates the number 


of cycles to failure. 


Residual stress Relaxation Exponents. 
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The name of a routine calling for integer-to-real conversion. | 
A variable subseript. 


‘The identi fying number of a block the blocks bali: 


numbered consecutively from r to NBLOCK. ° a * 


The name of a routine calling for real-to-integer conversion. 


: The nunber of steps daput., to the range patr counting 
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subroutine. 


Value controlling the WRITE statements. ao 


A counter. = 


7 Value. controlling entry ace the | ‘range par counting ne 
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An: index variable. . 


‘An index variable. 
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An index variable. 
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The ‘Adenti fying step number, tha steps being numbered. 
from 1 to NLEVEL. 7 = : : i ee 
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from 1 to NTYPE. 7 ae | 2 


A variable subscript. . 


“Value of +1 or 0, “as assigned for branch determination, 


Value of -1 or 0, as assigned for branch determination, 


An index variable. 


‘An index variable.. 


An index variable. 


Subvoutine..” 
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An index variable. 
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The total number of different types. 
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"Maximum applied stress. 


Minimum applied stress. 


Summation of damage for a ‘light. - 


(Summation of ENN/CYCLES). . 


Summation of R(N) , N=4,7. 


Summation of nR(N), N=4,7. n=4. 


Summation of wad N=4,7. 


Summation of ~ ‘ Ne4, 7. n=4, 
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Identification of the source of the SN data. 
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FATIGUE LIFE INPUT DATA FLR SEVERAL MATERIALS 
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31.7°° 
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